CO 

CD 

05 




NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 1963 



00 

CO 

CO 

CO 


! 

/ 


STRESSES IN AND GENERAL INSTABILITY OF 
MONOCOQUE CYLINDERS WITH CUTOUTS 
Vm - CALCULATION OF THE BUCKLING LOAD OF CYLINDERS WITH 
LONG SYMMETRIC CUTOUT SUBJECTED TO PURE BENDING 
By N. J. Hoff, Bruno A. Boley, and Merven W. Mandel 
Polytechnic Institute of Brooklyn 



Washington 
October 1949 



TECH LIBRARY KAFB, NM 




TECH LtPRARY KAFB, NFSfl 



00L5322 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 1963 


STRESSES IN AND GENERAL INSTABILITY OF 
MCNOCOQUE CYLINDERS WITH CUTOUTS 
VTH - CALCULATION OF THE BUCKLING LOAD OF CYLINDERS WITH 
LONG SYMMETRIC CUTOUT SUBJECTED TO PURE BENDING 
By N. J. Hoff, Bruno A. Boley, and Merven W. Mandel 

SUMMARY 


A strain— energy theory similar to one developed earlier at PIBAL 
(Polytechnic Institute of Brooklyn Aeronautical Laboratories) was 
established for the calculation of the buckling load in pure bending 
of reinforced mono cogue cylinders which have a symmetric cutout on the 
compression side and buckle according to general— instability patterns. 
The difference between the present and the earlier theories is the use 
of the axial wave length as an additional parameter whose value was 
determined from a minimum condition. The theory was applied to four 
cylinders which were tested earlier at PIBAL. Fair agreement was found 
between theory and the results of the experiments. 


INTRODUCTION 


Reinforced monocoque cylinders are subject to failure by a 
simultaneous buckling of frames, stringers, and sheet covering. This 
type of failure is known as general instability. The problem of the 
general instability of reinforced circular cylinders subjected to pure 
bending has been investigated in some detail by various authors (refer- 
ences 1 to 13 ) • These papers have dealt with complete cylinders, 
cylinders having symmetric cutouts on the compression side, and 
cylinders having side cutouts. 

Reference 12 presents a' strain-energy theory and the results of 
calculations for the buckling load of cylinders having a short 
symmetric cutout on the compression side. In the theory it was assumed 
that the general— instability bulge had the same wave length as the 
cutout . 

In this paper the general— instability buckling load of cylinders 
having long symmetric cutouts on the compression side is calculated by 
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strain— energy methods . The axial wave length of the "bulge is used as a 
parameter whose magnitude is found from the requirement that the "buck- 
ling load "be a mi nim um. Thus it is possible to deter min e whether the 
axl al wave length is the same as the length of the cutout. 

Test results showed that for some cylinders the bulge appeared to 
be antisymmetric with respect to the vertical plane of symmetry of the 
cylinder perpendicular to the axis. Therefore the calculations were 
carried out for the two different assumptions of deflected shape at 
buckling, one longitudinally symmetric, the other antisymmetric. 

In both assumptions for the shape of the bulge, the displacements 
at buckling were represented by the first seven terms of a Fourier 
series in the circumferential direction and by a few terms of a 
trigonometric series in the axi al. direction. In the case of longitu- 
dinal symmetry, the trigonometric series in x consisted of sine power 
terms and contained a parameter which had to be determined from the 
minimum buckling load condition. For longitudinal antisymmetry the 
series was a Fourier sine series in x which also contained a parameter. 
Boundary conditions were used to determine four of the seven coeffi— . 
cients of the Fourier series in <p for the circumferential variation 
of the shape while one remained undetermined as in all buckling problems. 
The remaining two coefficients as well as the circumferential wave 
length were calculated from the minimum buckling ■ load requirement. 

The strain energies considered were those due to bending of the 
rings in their own plane, radial and tangential bending as well as 
torsion of the stringers, and shear in the sheet. The extensional 
strain energy of the sheet was accounted for by adding an effective 
width of sheet to the rings and stringers. The external work was 
calculated on the assumption of a linear force distribution in the 
stringers, which is in better agreement with test results (reference 9) 
than a linear stress or strain distribution. 

The buckling load was calculated from the requirement that the 
strain energy corresponding to the transition from the unbuckled into 
the buckled shape be equal to the work done by the applied loads. The 
minimum value of the buckling load was found by assuming the circum- 
ferential wave length equal to some integral number of stringer fields, 
the axial wave length equal t<? some integral number of ring fields, and 
a numerical value for the parameter in the axial variation of the shape. 
The values of the two undetermined Fourier coefficients were calculated 
so as to make the buckling load a minimum, and with these two coeffi- 
cients fixed the value of the axial parameter was calculated in a 
similar mann er. The process was repeated until the three parameters 
reached steady values and the buckling load became a minimum. This 
value of the buckling load was compared with values obtained from other 
assumptions of axial and circumferential wave lengths and the lowest 
load was taken as the absolute minimum buckling load of the cylinder. 
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SYMBOLS 


a,a Q ,a^ 3 a 2 ,a 2 Fourier coefficients 

AjAeff cross-sectional area of stringer plus effective width 

of sheet 

b_,b 1 ,b 2 jb^ Fourier coefficients 


c 

C 

d 

e 

E 

®tan 

®red 

G 

Go 


G eff 

h 

i • 
Ir 


^StTy 


parameter 

geometric coefficient in torsional rigidity GC 

width of panel measured along circumference 

parameter 

Young’s modulus 

tangent modulus 

reduced modulus 
shear modulus 

shear modulus of sheet covering at zero compressive 
load 

effective shear modulus 
width of stringer 

index indicating position along circumference 

moment of inertia of ring section and its effective 
width of sheet for bending in its own plane 

moment of inertia of stringer section and its effective 
width of sheet for bending in the radial direction 
(about a tangential axis) 



b 


MCA W 1963 



3 


moment of inertia of stringer section and its effective 
width, of sheet for tending in the tangential direction 
(about a radial axis) 

index Indicating position along axial direction 


L 


length of bulge in .axl al direction 


distance between adjacent rings 


m number of rings involved in bulge 

m + 1 number of ring fields involved in bulge 


M 

Mcr 


n 


applied bending moment; function of n, a, and bj 
appearing in strain energy of bending in rings 

applied bending moment at buckling 

parameter defining length of bulge in circumferential 
direction 


Pl->P2 

p cr 

P i 

Q 

r 

R 


s 

S 

t 

U 

Ur 


polynomial functions of a and b 

force in one of the stringers at edge of cutout at 
buckling 

force in ith stringer 

function of x and c or e appearing in shear strain 
energy 

radius of cylinder 

function of n, <p, a, and b appearing in shear 
strain energy 

number of stringers involved in one— half of bulge 
total number of stringers in cylinder 
thickness of sheet covering 
strain energy 

strain energy stored in rings because of bending of 
rings in their own plane ; 
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^str 

r 

U str t 


*t 

Ush 

2w 

w n 

w r 

v t 


¥ 

x 

2a 


r 

6 


e 


cr 


9 


®r 


t 


0 

r 



strain energy stored in stringers "because of "bending 
about a tangential axis 

strain energy stored in stringers "because of "bending 
about a radial axis 

strain energy stored in stringers "because of torsion 
strain energy stored in sheet covering "because of shear 
effective width of sheet 
rotation of tangent to ring 

radial displacement of a point on a ring or a stringer 

tangential displacement of a point on a ring or a 
stringer 

work done "by external forces 

coordinate measuring distance along axis of cylinder 
cutout angle 

coefficients used in calculation of shear strain in a 
panel due to rotations and displacements of its 
comers 


shear strain 

distance of neutral axis from horizontal diameter of 
cylinder 

normal strain in a stringer 

"buckling strain of a panel of sheet covering 


angular coordinate measuring circumferential distance 
from edge of cutout 


function of n, cp, a, and "b appearing in expression 


for w, 


n 


function of n ; cp, a, and "b appearing in expression 
for w r 


function of 
for w. 


n^ 9^ a, and 


appearing in expression 
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TEE TIF,tfl',f( Tro SHAPE 


The shape of the "bulge at buckling is determined mainly by the 
radial deflections. It was observed in experiments conducted on 
cylinders designed to fail by general instability that, while most of 
the cylinders failed in a longitudinally symmetric bulge, there were 
some specimens which failed in what appeared to be a longitudinally 
antisymmetric bulge pattern (reference 14). For that reason, these 
two types of deflected shape were considered. Typical deflection 
patterns of stringers and a ring are shown in figure 1. 


SIMMETRIC DEFLECTION PATTERN 


The expression chosen to represent the radial deflections that 
are symmetric longitud in al l y is: 




= sin^(itx/L) + c sin^(jtx/L) & 


( 1 ) 


where 


$ r = /'aQ + a-^ cos nrp + a£ cos 2ncp + a^ 


cos 3nrp + 


sin nip 




(la) 


provided 


0 ^ cp — jr/n 
0 < x < L 


Also, 


w r = 0 


for cp > jr/n and/or x < 0, x > L. The notation and sign convention 
are shown in figure 2. 

The deformations of the rings were assumed to be inextensional. 
The condition for inextensional deformations is: 


w r = — Sw-fc/ckp 


( 2 ) 
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Therefore the tangential deformations are given by: 

w t = [sin 2 (rtz/L) + c sin 6 (jtx/L)]$t (3) 

where 

= — (l/n)|aocp + aq sin mp + sin 2ncp + ^3/3) sin 3ncp - 

bq cos ncp — (^2/2) cos 2ncp — (b^/3) cos 3ncpJ (3a) 


provided 

0 < cp < it/n 
0 < x < L 


Also, 


w t = 0 

when cp > jr/n and/or x < 0, x > L. The arbitrary function of x 
which would normally appear in equation (3) as a result of the inte- 
gration of equation (2) is zero if the displacements are symmetric 
about a vertical diameter. 

If it is required that there be a smooth transition between the 
bulge and the undistorted part of the cylinder at qp = n/n, then; 

(1) The tangential displacements must vanish: 

w t = ° (4a) 

when (p = jt/n for all values of x 

e> 

(2) The radial displacements must van i sh: 

w r = 0 (4b) 

when cp = Jt/n for all values of x 

(3) There must be no sudden change in the direction of the tangent: 

dw r /&P = 0 (4c) 

when cp = ir/n for all values of x 
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(4) There must "be no sudden change in the curvature: 


8 2 w r ^3cp 2 = 0 


(4a.) 


when qp = n/n for all values of x. Also, in order to have a smooth 
transition between the "bulge and the undistorted part of the cylinder 
at x = 0 and/or x = L: 

(5) The tangential displacements must vanish: 

w t = 0 ’ (he) 

when x = 0 or 1 for all values of cp 

( 6 ) The radial displacements must vanish: 

v r = 0 (4f ) 


when x = 0 or L for a ll values of (p 

(7) There must "be no sudden change in the direction of the tangent: 


dw r fdx = 0 


(*«) 


when x = 0 • or L for all, values of cp. Conditions ( 5 ), ( 6 ), and (7) 
are automatically satisfied hy the expressions in the assumed deflected 
shape for all values cf the parameter c. The remaining four conditions 
establish four relationships between the Fourier coefficients and make 
it possible to determine any four of the coefficients in terms of the 
remaining three. If a^, a^, and b^ are retained as the basic param- 
eters, and the notation 


a l/ a o = a 
bi/a Q = b 


(5) 


is used, then equations (la) and ( 3 a) become: 


w r = a Q [j 3 in 2 (jtx/L) + c sin 6 (jtx/L)] [l + a cos ncp + (1.6a — 1.8) cos 2ncp + 
( 0 . 6 a — 0 . 8 ) cos 3 ncp + b sin ncp + ( 3 . 2 b + 3 . 6 n) sin 2 ncp + 

(l. 8 b + 2.4ir) sin 


( 6 ) 
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w t = -< a o/i 1 )[sin 2 (jtz/L) + c sln^(«x/L)] [ncp + a sin nrp + 

(0.8a — O. 9 ) sin 2nrp + (0.2a — 0.2 666 . . .) sin 3 nc P — 
b cos np — (1.6b + 1.8 jt) cos 2np — (0.6b + 0.8ir) cos 3rapl ( 7 ) 


CALCULATION OF STRAIN ENERGY FOR THE 
SYMMETRIC DEFLECTED SHAPE 
Strain Energy Stored in Rings 


Tbe strain energy stored in one— Half of any one ring is; 


U 


" ( EI r/ 2r3 ) f [ w r + ( /8cp 2 )] ‘ 


Jo 


dp 


( 8 ) 


If the value of w r is substituted from equation (6) and the strain 

energy is summed up over all the rings, the following expression is 
obtained: 


U r = C 1 / 2 ) 2Z ( E I r /r 3 ) [sin^ ( it j /L ) + 

J=1 


c sin' 


’(rtj/L )] 1 



+ 



(9) 


where m is the total number of rings included in- the wave length L. 
The integration yields a result in closed form. If al 1 the rings have 
the same bending rigidity EI r , the total strain energy for one— half 
of the cylinder is : 


u r = ( a o 2 / 2r3 )( EI r M ) YZ + c sin 6 (jtj/L)J 

J=1 


( 10 ) 



10 


HA.CA TN 1963 


where 

nM = [rt + 10.053096(l - 9n 2 ) + 206.01005(l - 4n 2 ) 2 + 

90.303387(1 - 9n 2 ) 2 - l8.095573(l - 4n 2 )(l - 9n 2 )] + 
a|— 9.0477868(l - 4n 2 ) 2 - 1.5079645(1 - 9n 2 ) 2 + 

30.159289(1 - n 2 )(l - 4a 2 ) + 18.095573(1 - 4n 2 )(l - 9a 2 )] + 
b [4(1 - a 2 ) + 2.4(1 - 9n 2 ) + 113.69784(l - 4a 2 ) 2 + 

42.636690(1 - 9a 2 ) 2 + 2.4(1 - a 2 )(l - 4a 2 ) - 
3.68(1 - 4a2)(l - 9n 2 )] + a2[(jt/2)(l.- a2) 2 + 4.0212386(1 - 4a2) 2 
0.5654867(1 - 9a 2 ) 2 ] + b 2 [(*/2)(l - n 2 ) 2 + l6.084954(l - 4a 2 ) 2 + 
5.08938(1 — 9n 2 ) 2 ] + ah j^6.4(l — a 2 )(l — 4a 2 ) + 

3.84(l - 4a 2 ) (l - 9n 2 )] (11) 

It is possible hut cumbersome to sum up iu closed form the trigonometric 
fuactioas describiug the deflected shape ia the axial directioa as 
given in equation (10). It was found more convenient to carry out the 
summation numerically. 

Strain Energy Stored ia Stringers 

The strain energy stored in the stringers because of heading ia 
the radial direction is: 

^str r = (l/ 2 ) y I'lstrp [ (8 2 w r /8x2) dx (12) 

1 Jo 

where the summation is extended over all the stringers involved in the 
bulge in one-half the cylinder. Substitution of the value of w r from 
equation (6) into equation (12) and integration gives: 
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P L 


^str = C 1 / 2 ) H EI str ($ r ) 2 ( rt /L) if - 

r i r 


[2 cos 2 (itx/L) - 2 sin 2 (jtx/L) + 


0 


30c sin^mx/L) cos 2 (jtx/l) — 6c sin^(icz/L)J 2 dx 
u str = (^/L 3 )^ + (15/8)0 + (I^l/I28)c 2 ] EI str ($ r ) 2 


(13) 


The strain energy stored in the stringers because of bending in 
the tangential direction is; 


U str t - t 1 / 2 ) H EI strt 


(^t/ax 2 )^ 


dx 


(ik) 


0 


where the summation is extended over all the stringers involved in the 
bulge in one— half of the cylinder. With the aid of equation (7) this 
strain energy becomes: 

u str t = (" 4 / l3 )[ 1 + (15/8)o + (Ma/I28)c 2 ] JZ EI strt («f t ) 2 (15) 
The strain energy stored in the stringers because of torsion is: 

n t = (1/2) JZ 00 

i 

and again the summation is extended over all the stringers involved in 
the bulge in one-half of the cylinder. In this equation (l/r) (d 2 w r /8xckp) 
is the unit angle of twist for the stringers. In the expression for the 
Saint Tenant torsional rigidity 

C = O.lW 1 ' (16a) 

since the test specimens were provided with square section stringers of 
edge length h. Differentiation gives : 


(l/r) 2 (8 2 w r /dx8cp) 2 dx 


(16) 


d^p/dxckp = (8/dx) [sin 2 (jtx/L) + c sin^(nx 


(IT) 
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where 


4 


r 


t 


a 0 n£— a sin ncp — ( 3 . 2 a — 3 . 6 )sin 2 nrp — ( 1 . 8 a — 2 . 4 )sin 3 ncp + 


h cos ncp + (6.4b + 7 . 2 jt)cos 2 ncp + (5.4b + 7 » 2 jt)cos 3 nrp] ( 17 a) 
This gives for the strain energy of torsion: 

U t = (fl^GC/Lr^ ) [( l/4 ) + (15/32 )c + (189/512)c. 2 ] II ($ r ’)- 2 (18) 

where the summation is over one-half the cylinder, as before. Since the 
variation of the torsional rigidity caused "by the different amounts of 
effective width of sheet is small when calculated according to the 
Saint Tenant theoiy, it was considered permissible to assume GC a 
constant. 

The values of 4 r , 4^, And 4 r ’ vary from stringer to stringer. 
Also I s trr ^strt Yar J around the circumference because the 

effective width to be added to each stringer changes. Therefore the 
summations indicated in equations ( 13 ), ( 15 ), and ( 18 ) had to be 
evaluated numerically for each cylinder investigated. 


Strain Energy of Shear Stared in Sheet 

The shear strain energy per unit volume was taken as one— half the 
average effective shear modulus multiplied by the square of the average 
•shear strain 7 for each panel of sheet. The average shear strain '7 
was calculated from the relative displacements of the corners of the 
panel, as shown in figure 3* Then the total strain energy of shear 
stored in the sheet is: 


U sh = (1/2) H 7 2 G eff Ljtd 


(19) 


where Lqtd is the volume of one panel and the summation extends over 
al 1 the panels involved in the bulge in one— half the cylinder. The' 
effective shear modulus G eff depends on the geometric and mechanical 

properties of and the average normal strain in the panel. Its value was 
taken from the empirical curves established earlier at PIBAL and 
presented in figure 24 of reference 15* 
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The average angle of shear 7 was calculated from the equation: 

7 - (arAx^.j - ”r 1+ltJ ~ "r ljJ+1 * * 

(N^i)Kj + Nrt.j - Xl+i ' + 

+ - %, J+ i - % + i jJ+ i) (20) 

where the subscript i refers to the circumferential location, and the 
subscript j, to the axial location of the comers of the panel. 

• The rotation w n of the tangent of the ring is given by the 
relation: 


w n = (l/r)(dw r /dcp) = (l/r)®! 


( 21 ) 


Tn reference 1 6 the values of the numerical factors a^, a^., and 
were determined. 


ctp = d/lQr = 2 rc/lCS 

= — (l/2)[l - 0.01666 . . . (d/r) 2 ] 
= — (r/d)(0.5 + c^) 


Substitutions yield: 

m 

u sh = (tdG Q /2Li) ZZQjZI (Geff/Go)i% 


s— 1 


( 22 ) 


( 23 ) 


where Q is a function of x only, and E, a function of cp only, 
given by the relations : 

Q = ^|sin 2 [jtj/(m + 1)] + c sin^[jtj/(m + 1)]|- |sin 2 [n(j + l)/(m + 1)]] 
c sin^[jt(j + l)/(m + l)]j^ 


( 24 ) 
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E = -*r 1+1 )+ KIK + *t 1+1 ) + 

(l°* | ^X»r ± ' + »r 1+1 ’)] 2 (25) 


It is possible to get a result in closed form for the summations Q 
and E "but for convenience the summation was performed numerically. 


WORK DONE BY EXTERNAL FORCES * 


It was observed in the experiments described in reference 9 that 
the stress distribution was not linear in the cutout portion of the 
cylinder, although the deviations from linearity were not large as a 
rule. A good approximation to the experimental curves was obtained by 
assuming a linear force distribution, which is not equivalent to a 
linear stress distribution because of the varying amounts of effective 
width of sheet added to the stringer sections. A comparison of the 
strain distribution calculated on the assumption of linear force distri- 
bution with the experimental strains is given in figure 4 of reference 12. 
The expression used for the calculation of the force acting on the ith 
stringer is: 


P-, = P. 


cr 


[5/] 


r + cos 


(a + 2*i/s)]/(5/r + cos a) 


(26) 


where P„^ is the compressive force acting upon the stringer at the 
cr 

edge of the cutout, and 5 . is the distance of the neutral axis from the 
horizontal diameter of the cylinder. 

Equal and opposite forces are assumed to be acting at the x = 0 
and x = L ends of each stringer. The distance between the points of 
application of these forces shortens at buckling. The work done by the 
forces is equal to the summation of the forces times the shortening of 
all the stringers involved in the bulge in one— half the cylinder: 


¥ = 



(dw-fc/cix) 2 dx 


(27) 
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Substitutions and integration give; 

¥ = 1/4) + (15/32)c + (189/512) c 2 ]^^ 2 + « t 2 ) 

¥ = (* 2 /L)[(lA) + (15/32)c + 

(l 89 / 5 l 2 )c 2 ]P Gr E(*i/Par) (* r 2 + ® t 2 ) ( 28 ) 

The sum m ation in the right-hand member of equation (28) was carried out 
numerically. 


ANTISYMMETRIC DEFLECTION PATTERN 


The expression chosen to represent the radial deflections tha t are 
antisymmetric longitudinally Is: 

w r = |sin (itx/L) + e sin (h-icn/L) — 

[(1/3) + (2e/3)] sin (6nx/L)J& r * (29) 

for 

0 < cp < jt/n 
0 < x < L 


where <& r is given by equation (12). Also, 


w_. 


= 0 


for cp > Jt/n and/or x < 0, x > L. The expression for the tangential 
deflections, again on the assumption of inextensional deformations, c an 
be obtained from equations ( 2 ) and ( 29 ) and is given by 

w^ = -|sin- (2jtx/L) + e sin (4 «x/l) — 

[(1/3) + (2e/3)] sin ( 6 «x/L)|® t 


(30) 
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for 

0 < cp < it/n 
0 < z < L 

where is given by equation (3a). Also., 


w t = 0 


for qp > rt/n and/or x < 0, x > L. These expressions for w r and w + 

must satisfy the boundary conditions at <p = jt/n given in equations (4a) 
through (4g). 

It can be shown that the expressions for w r and w^ given by 
equations (29) and (30) satisfy the given boundary conditions for al 1 
values of the longitudinal parameter e, as well as the parameters a 
and b. In addition, it can be verified that equations (29) and (30) 
satisfy the requirement for antisymmetry in the axial direction, namely, 

r (x) r (L— x) (31) 

for any given value of cp. 

CALCULATION OF STRAIN' ENERGY FOR THE ANTISYMMETRIC 
REELECTED SHARE 


The strain energies are derived in the same way as those for the 
symmetric deflected shape. They are listed below with reference to the 
equations from which they are derived. 


Strain Energy Stored in Rings 

From- equations (8) and (29) the strain energy of bending at the 
rings involved in the bulge in one-half the cylinder is : . 

m 

^r = (ao 2 / 2r3 ) EI rM JZ 

[(1/3) + (2e/3)] sin (6ttj/L)J 2 (32) 

where M is given in equation (11). The summation is taken over the 
rings involved in the bulge and was evaluated numerically. 


■ 


sin (2jrj/L) + e sin (4irj/L) — 
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Strain Energy Stored in Stringers 

From equations (12) and (29) the strain energy due to radial 
bending of the stringers is given by 


U str r (*Vl 3 )(iK> + I44e + 208e 2 )E^EI strr ($ r ) 2 (33) 

From equations (14) and (30) the strain energy due to tangential 
bending of the stringers is 

U str t = (^A 3 )(^0 + Ikke + 20Qe 2 )^Zm Btr ^ t ) 2 (34) 

Equations (1 6) and (29) give for the strain energy due to torsion 
of the stringers: 

tJ t = + 4e + 8e 2 ) (GC/r 2 ) J ~ ($ r *) 2 (35) 

where C is given by equation (l6a). The summations for all stringer 
strai n energies were evaluated numerically. 


Strain Energy of Shear Stored in Sheet 

From equations (19), (21),, (29), and (30), the strain energy of 
shear in the sheet is; 

m a— 1 

U sh = ( tdG o/2l!) YL Qj ZZ (&eff/G-o)i% ' (36) 

J=0 i=0 

where Q is given by 

Q = ^£sin [ 2 *«j/( m + 1)3 + © sin [4rtj/( m + 1)] - 

[(1/3) + (2e/3)] sin [6*j/( m + 1)]J- | sin [2«( J + l)/(m + 1)] + 
e sin [4ir( j + l)/(m + 1)] - [(1/3) + 

(2e/3)] sin [6n(j + l)/(m + 1)] ^ 

and B is given by equation (25). The summations y~Q and y~ B were 
evaluated numerically, ‘ 
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WORK DONE BY EXTERNAL FORCES 

From equations ( 27 ), (29), and. (30) the work done By the external 
forces is: 


W = («2/L) P cr (2 + he + 8 e 2 )Z^(Pi/P cr )( 3 ) r 2 + <b*) (37) 


CALCULATION CF BUCKLING LOAD 


The Buckling condition is 


U r + ^str r + u str t + u t + U sh = W 


( 38 ) 


where the values of the strain energies and the work must Be taken from 
equations (10), (13), (15) j (18), (23 ) J and (28) for the longitudinally 
symmetric Buckling shape, or fro^i equations (32)^ (33), ( 3 ^), ( 35 ), (36), 
and ( 37 ) for the longitudinally antisymmetric Buckling shape. 

Equation (38) was solved for P cr contained in W in the following 
way. 


Integral numbers of s and (m + 1 ) were chosen for the number 
of stringer and ring fields included in the Bulge. For these values 
the M, and E functions reduced to quadratic expressions in a 

and B. Next P cr was assumed. This permitted the calculation of the 
effective widths of sheet acting with the stringers, the moments of 
inertia of the stringers, and the values of the effective shear modulus 
in the sheet panels. 

Finally, a trial value was chosen for the parameter c (for the 
symmetric deflected shape) or e (for the antisymmetric deflected 
shape) in the expressions for the- longitudinal variation of the 
Buckled shape. All necessary summations were then carried out. Substi- 
tution of al 1 the results in equation (38) made it possible to obtain a 
solution for P cr in the form; 


P cr = Pl( a / b )/?2( a *' b ) ( 39 ) 

where p^_ and P2 are' quadratic expressions in the parameters a 
and. B. Minimizing this expression for P Qr with respect to a and B 
is equivalent to setting 
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P cr = Pi( a *‘ b VM a * 1> ) = (Sp!/aa)/(ap 2 /Sa) = (dpi/db)/(dp 2 /db) (40) 


where the partial differential coefficients of p 1 and p 2 are linear 

functions of a and b. Equations (40) give three relations 
between P cr , a, and b, which may be solved by a trial— and— error 

procedure. With the aid of an assumed value for P cr , a and b were 
calculated from the linear equations . These values of a and b were 
then substituted into the original quadratic expressions in equa- 
tion (39), yielding a calculated value for P cr . This procedure was 
repeated with new assumptions for P cr until the value calculated was 
reasonably close to the one assumed. These values of a and t were 
then substituted in the original expressions for the strain energies, 
thus reducing the 0 M, and R functions to quadratic expressions 
in c or e, depending on which type of deflected shape was under 
consideration. The necessary summations were carried out again, and 
substitution of the results in equation ( 38 ) gave a solution for P 
in the form cr 

P cr = P 3 ( c )/P 4 -( c ) (4la) 

or 


p cr = P 5 (®)/P 6 ( e ) 


(41b) 


where p and p^ are quadratic expressions in the parameter c, 
and pc; J and pg are quadratic expressions in the parameter e. 

Minimizing this expression for P cr with respect to c or e is 
equivalent to setting 

P cr = P 3 ( c )/P4( c > = (dP 3 /dc)/(dp 4 /dc) (42a) 

or 

p cr = P 5 (e)/pg(e) = (Sp 5 / / 8 e)/( 8 p 6 /Se) (42h) 


where the partial differential coefficients of p^ and p^ are linear 
functions of c, and the partial differential coefficients of p^ 
and pg are linear functions of e. These equations can be solved by 
the same trial— and— error procedure that was used to solve equations (40) 
for a, b, and P cr * If the value of the parameter c (or e) was 
different from the trial value assumed at the very outset when the 
parameters a and b were determined, then the new value of c (or e) 
was substituted in the expressions for the longitudinal variation of the 
buckled shape and the minimization procedure of equation (40) was 
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repeated in order to get new values for a and 1. The process was 
continued until the values of all parameters reached reasonably steady 
values and the value of P cr approached a minimum . 

This procedure involving separate minimizations for the longitudinal 
and circumferential parameters was used because minimization with respect 
to all the three .parameters at the same time would have yielded nonlinear 
equations in the parameters. It was considered more convenient to repeat 
the minimizations than to solve the nonlinear equations. 

The procedure was carried out for different numbers of ring and 
stringer fields included in the bulge. The smallest buckling load 
obtained in this way was considered the true buckling load. 

When the value of P cr obtained in these calculations differed 
materially from that assumed at the outset* the moments of inertia and 
effective shear modulus had to be calculated again and the entire 
procedure repeated. 

Details of the procedure may be seen from the numerical example 
given in the appendix. 


COMPARISON OF THEORY AND EXPERIMENT 


Numerical calculations were carried out to obtain the buckling 
load* for the axially symmetric deflected shape* of the cylinders shown 
in figure k. The buckling load for the axially antisymmetric deflected 
shape was determined only for cylinder 75* Typical deflected shapes of 
a ring and stringers are shown in figure 1 and the results of the calcu- 
lations are given in table 1. The experimental bending moments at 
buckling were taken from reference llj-. 

The theoretical buckling loads were consistently higher than the 
experimental values* and the deviations between theory and experiment 
seemed to increase systematically with increase in ring size. Compar- 
ison of the theoretical buckling loads for axially symmetric and anti- 
symmetric deflected shapes shows that the symmetric shape gave a lower 
value but the difference was small. This indicates that the assumption 
of a symmetric shape was closer to the true deflected shape than the 
assumption of antisymmetry. 

The axial wave length predicted by theory was found to be different 
from the length of the cutout. 

Comparison of results obtained with the Jj-5° cutout cylinders of 
reference 12 with cylinder 72 of this series which has rings of the 
smallest size shows a decrease in the 'deviation between theory and 
experiment from an average value of 33-9 percent to 18.1 percent. It 
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appears that the assumption of a more complex deflected shape incorpo- 
rating- more free parameters was the cause of this improvement. However, 
in the case of the heaviest rings, even the present assumptions are not 
sufficiently general to give a "buckling load in reasonable agreement 
with experiment. 


.CONCLUSIONS 


A strain— energy theory was developed for the general— instability 
buckling load of reinforced circular monocoque cylinders having a long 
symmetric cutout on the compression side and subjected to pure bending. 
When the theory was applied to four test cylinders of the experimental 
reports, the buckling loads obtained were 18 . 1 , 26. 5 , 25 . 6 , and 45-1 per- 
cent higher than the experimental values. The average deviation was 
28.8 percent. 


Polytechnic Institute of Brooklyn 
Brooklyn, N. Y., July 12, 1948 
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APPENDIX 

NUMERICAL EXAMPLE 


In this appendix details of the calculations performed in deter- 
mining the "buckling load for PIBAL cylinder 73 are shown for the case 
of a longitudinally symmetric deflected shape. The geometric and 
mechanical properties for this cylinder are: 


Sheet thickness, in 0.012 

Radius, r, in 10 

Distance between adjacent rings, L^, in 3*856 

Stringer spacing along circumference, d, in 3*927 

Number of stringers in full portion of cylinder, S 1 6 

Angle of cutout, 2a, deg . . • . . . ^5 

Stringer cross section, in 3/8 by 3/8 

Area of stringer cross section. A, sq in 0.140625 

Torsional rigidity, C, in. 1 * - 27.686 X 10 ~^ 

Ring cross section, in l/4 by l/P 

Moment of inertia of ring cross section, I r , in.^ . . . 7^9 . 36k X 10~° 

Modulus of elasticity, E, psi 10. 5 x 10° 

Shear modulus, G = G 0 , psi. 3.9 x 10^ 


At the outset, the assumed value of P cr was taken as 53OO popnds. 
The shift of the neutral axis was calculated to be 1.14-52 Inches on the 
basis of a linear force distribution for pure bending. This permits the 
setting up of the following table which gives the effective width 2w 
of sheet to be added to the stringers, the moments of inertia of the 
stringer cross sections, the effective shear modulus G e ff, and the 
reduced moduli of elasticity. 



[s = 1.1452; € or - 3.3 x lcr+j 






( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 

( 9 ) 

( 10 ) 

(U) 

( 12 ) 

(13) 

( 14 ) 

B 


€ A 

G 



A eff 

2 w 

I atr r 


■^tan 

*r«d 

B 

e 

e /*or 

0 offA) 

0 

5300 

5.048 x 10~4 

33.5 X icr 4 

0.1506 

1.281 

19.68 X lcH* 

25.O x 10“^ 

9,79 x 10 6 

H 

O 

£ 

X 

I-* 

0 + 1/2 

29.3 x lcr^ 

8.88 

0.475 


4198-9 

3-999 

25.3 

.1582 

1.465 

21.76 

48.0 

10.5 

10.5 

1 + 1/2 

20.4 

6.18 

.497 

2 

2550.9 

2.429 

14.8 

■ 1637 

1.923 

23.32 

87.0 

IO.5 

10.5 

2 + 1/2 

9.0 

2.73 

.620 

3 

606.9 

-578 , 

3-1 

.1878 

3.927 

23.92 

623.0 

10.5 

10.5 

3 + 1/2 


— 

1.000 


In this table., column (1) refers to the stringer station. Column (2) is obtained from, 
equation (2 6), column (3), by dividing oolumn (2) by K “ 10 . 5 X 10°psi. Columns (4), ( 5 ), 
and ( 6 ) oan be most conveniently filled In through the use of a previously drawn curve of . 
the strain g against tha arsa - ^- q ±x of* stringer and. ©iTfactivs width, o omi) inat i on • (A 
curve of this type was used in the present calculations, and was constructed with the aid of 
the following fo rmu la for the effective width: 

(Al) 

derived in referenoe 1 . ) 


2 v =» (l/e)(d/r) 


0.3t + 1-535 [(t/d)(er - O^tlr 1 / 2 ] 2 / 3 } 


J 


Columns (7) and (8) give the moments of inertia of the stringers plus their effective 
width of a curved sheet calculated from the equations : 

Istr r = (hVl2) + - fa + t V 2 - [(2w) £ /2l|r]) + 0 .8 [(2v) 2 /24r] 2 2wt 

(l/2vb) + (1/A) 1 J 


(A2) 


Istr't = ( h Vl2) + (1/I2)(2w)3t 


(A3) 


ro 

U) 
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where h is the width of the stringer. These equations were also 
derived in reference 1. 

If the strain in the stringer is above the proportional limit, it 
is necessary to use the von Earman reduced modulus of elasticity for 
calculating the bending rigidity of the stringers. Columns (9) and ( 10 ) 
give the values of the reduced moduli based on the curves of 
reference 17 • 

Column ( 12 ) gives the value of the strain in the middle of each 
panel from which the effective shear modulus can be determined. On the 
basis of a buckling strain of the panels e cr = 3.3 X 10 -^, column (13) 
gives the ratio e/e cr for each panel. The values of Cr 0 ff/Cr o in 
column (lk) are obtained from figure 2 k of reference 15. 

For the evaluation of the 4 >^, M, and E functions, the number s 
of stringer fields included in the bulge in one— half the cylinder was 
assumed to be four. The corresponding value at n can be obtained 
from equation (Al)t 

n = S/ 2 s (Ak) 

Only integral numbers of stringer fields were used in the summations, 
so that the angle cp was replaced by its equivalent (2jri/s), 
where ( 2 n/S ) is the angle subtended by one stringer field, and i = 1, 
2, 3, . . . denotes the circumferential location of each stringer. 

The trigonometric functions appearing in the expressions for <D r , 
and ® r f are given in the following table, together with the coeffi- 
cients of the trigonometric functions. 


O = k; n = 2] 


i 

Constant 

“P =IT 

_ iti 

cos 

s 

_ 2 jt± 

cos 

s 

cos 35 * 
s 

sin — 

S 

, 2 jti 

sin 

8 

sin^ 

8 

0 

1 

0 

1 

1 

1 

6 

0 

0 

1 

1 

0.785398 

0.707107 

0 

-0.707107 

0.707107 

1 

0.707107 

2 

1 

1.570796 

0 

-1 

0 

1 

0 

-1 

. 3 

1 

2.35619k 

-0.707107 

0 

0 

• 

3 

H 

O 

0.707107 

-1 

0.707107 

Multiplier 

1 

0 


1 . 6 a 

0.6a 

t 

3 . 2 b 

1.8b 

for ® r 




- 1.8 

— 0 . 3 a 


3-6 . 

2 .k 

Multiplier 

for 

0 

-0.5 

0.5b 

0 0 
Coo? 

0.3b 

O.k 

-0.5a 

— O.ka 
0 .k 5 

—0.1a 
O.I33 • • • 

Multiplier 
for Op’ 

0 

0 

2 . 0 b 

12 . 8 b 

lk.k 

10 . 8 b 

lk.k 

- 2 . 0 a 

— 6 .ka 
7.2 

— 3 - 6 a 

k .8 







































The polynomial <t r for 1=0, for example. Is obtained by multiplying, column by 
column, the expressions in the row for 1 = 0 by the mul tipliers for 4 r and adding the 
resulting products. The results for $» r , 4^, and 4 r f are given in the following table. 


1 

®r 

H 

♦ 1 

0 

-1.6 

3-2a 

0 

4.0840705 

0 

1.6b 

90.477868 

0 

25.6b 

1 

18.2068803 

0.2828428a 

5.1798996b 

-0.7369949 

-0.464264a 

0. 1434214b 

-21.3946536 

-10.3598a 

-6.2225416b 

2 

-4.7398224 

' 

—1.6a 

-0.8b 

-3.7461647 

—0.4a 

-0.8b 

-50.0389342 

1.6a 

-12.8b 

3 

-5.543958 

-0.2828428a 

-1.2201004b 

-0.6452392 

-0.0242642a 

-0. 1414214b 

28.1828808 

2.4402008a 

6.2225416b 


The value of M oan be calculated with the value of n = 2 corresponding to 8=4 
and is 


M = 73562*5088 + 3487.1679a + 37945.9796b + 805.8186a 2 + 4933.8711b 2 + 1151. OOOab (A5) 

The functions (4 r ^ ^t^ + ^t-j^)* anii ^r^* + n^eied- for the calculation 

of B can be determined by adding or subtracting the polynomials In adjaoent rows for each 
column of the above table. However, It must be remembered that ai 1 4 functions are zero 

at i = 4. The results are not given here. For the values of r =■ 10 inches 
and d - 3.927 inohes equations (22) give: 
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a* = 0.03927 
a t = -0.49871 > 

ct^ = -O.OO3272 


(AS) 


Substitution in equation (25) gives the values of B which are listed 
in the following table.. 


[a r = 0 . 03927 ; a t = - 0 . 49871 ; a n = -O.OO3272] 


i 

E 

0 

-O.80694216 

0 . 2842619 a 

-0.6423940b 

1 

I.2486636 

0.5269357a 

0.0700127b 

2 

2.1360358 

0.26904896a 

0.4457607b 

3 

O.5765488 

0 . 0263752 a 

0.1266206b 


The $ 2 , Mj and E functions appear in the strain— energy expres- 
sions in quadratic form, but the squares of the polynomials are not 
given here. 

The next step in the calculations is to assume a value of (m + 1 ), 
the number of ring fields involved in the failure . The value chosen 
was (m + 1) =10 and the summations appearing in equations (10) 
and (23) were evaluated numerically on that basis. These results as 
well as the quadratic functions of c appearing in the expressions- 
for the strain energy of the stringers and the external work are given 
in the following table. 


. [(m + 1) = 10] 


Factor for strain 
energy of — 




Eing bending 
Stringer bending 
Stringer torsion 
Sheet shear 
External work 

3.750 

1.0 

O.25 

0.4774575 

O.25 

5 . 468750 c • 

1.875c 

0 . 46875 c 

0.8952328c 

0 . 46875 c 

2.2558595c 2 

3 . 4453125 c 2 

0 . 3691406 c 2 

0 . 6753475 c 2 

0 . 369 l 406 c 2 
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Finally a trial -value of the parameter c was chosen. In this 
numerical example the value of c used is not the ori ginal trial value 
"but one which was obtained after several trials. Substituting the * 
value c = -1.477 in the quadratic functions of c makes it possible 
to evaluate the various strain energies in the form of quadratic 
expressions in a and b for each value of i. 

The quadratic expressions of a and b in . <t 2 , M, and R are 
multiplied by the following factors from equations (10) , (13), (15 ), 
(l8), (23), and (28 ) , which have been multiplied by L /a 0 2 n 2 f Qr 
convenience of calculation: 

For ring bending: 

^EI r L/2Ttr 3 )(3.750 + 5-^6875c + 2.2558595c 2 ) = 9.IO99O (A7) 

For stringer bending: 

(e^/lXi.O + 1.875c + 3-^53125c 2 )(l stri ) = 400j399.8l stri (a 8) 

For stringer torsion; 

(Gc/r 2 )(o.25 + 0.46875c + 0.3691406c 2 ) = 39.16829 (A9) 

For shear of sheet: 

1 [o 0 tr(m + l)/ns](o eff /G 0 ) i (o.477 1 <-575 + 0.8952328c + 0.6753475c 2 ) 

= 58 J 482 . 6^G e f f /Go) i (A10) 

For external work: 

( p i/ p cr)(°* 2 5 + 0.46875c + 0.3691406c 2 ) = 0.362752^^) (All) 

Multiplication of the $ 2 , M, and R functions by these multipliers 
and addition of the results make it possible to solve for P cr : 

^ _ 1,804,446 + 147,812a + 935,^31b + 36,451a 2 + 23,960ab + 130,546b 2 

110.0681 + 2.74494a + 61.82786b + 4.27776a 2 + 8.93153b 2 + 1.39202ab 


(A12) 
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Differentiation with, respect to a and 1 , respectively, yields: 


> _ 1^7,812 + 23,960b + 72,902a 

cr 2 . 71^9!*. + 1.39202h + 8.55552a 

■ 

, 935^31 + 23,960a + 261.0921 

cr 61.82786 + 1.39202a + 17.863061 


(A13) 


These equations were reduced to two linear equations ly assuming a 
value of P cr = 5500 and clearing fractions. The two linear equations 
were solved simultaneously for a and 1: 


a = -3.019 

= - 3 - 35 ^ 


(Aik) 


Bulstitution of these Values in equation (A 12 ) gave a value 

of P cr = 5507 pounds. This result was considered sufficiently close 

to the value assumed originally. 


The next step was to substitute these numerical values of a 
and 1 into the $ M, and E functions in order to allow a mini- 
mization with respect to the parameter c. However, the <S M, and E 
functions happened to have leen previously evaluated with a = —3-07 
and 1 = — 3.36. Since small changes of the parameters have little 
effect on the buckling load, and since these values are very close to 
those given in equations (Al 4 ), the latter values were used. 


The numerically evaluated <t>^, M, and E functions were then 
multiplied ly the remaining factors for the various strain energies 
from equations (10), (13), (15), (l8), (23), and (28), exclusive of the 
quadratic expressions in c. This gave as multipliers for the c 
functions : 


For ring lending: 

l6l,8ll.8 

For stringer lending: 

19 , 103.99 

For stringer torsion: 

123 , 458.2 

For shear of sheet: 

29,527.52 

For external work: 

136.2066 


(A 15 ) 
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Multiplication of .the quadratic expressions in c by these factors and. 
addition of the results made it possible to solve for P cr as a 
function of c : 


671,407 + 1,006,097c j- 498,239c 2 
34.0517 + 63.8468c + 50.2794c 2 


(JUS) 


Differentiation with respect to c as in equation (42a) gives: 


. = 1,006,057 + 996,478c 

cr 63.8468 + IOO.5588C 


(A17) 


Assumption of P = 5440 yields a value of c = —1.4757 from 
equation (A17)« Substitution of this value of c in equation (Al6) 
gives P cr = 5446. Since this value of c is approximately equal to 
the trial value c = —1.477 assumed at the beginning, the value 
of P cr obtained is considered the minimum for the number of ring and 
stringer fields (s =4, (m + 1) = 10) under consideration. 


The entire procedure was repeated for other integral numbers of 
ring and stringer fields included in the bulge, and the lowest value 
of P cr was taken as the minimum value, and hence the value corre- 
sponding to the buckling load, for the cylinder. 


The final results of the procedure are listed in table 1 for the 
four cylinders investigated. The cylinders are shown in figure 4. 
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Figure 3.- Notation and sign convention for shear in a panel. 
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Figure 4.- Monocoque test cylinders. 
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